Evidence of a bond-nematic phase in LiCuVC>4 
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Polarized and unpolarized neutron scattering experiments on the frustrated ferromagnetic spin- 
1/2 chain LiCuVCU show that the phase transition at .Hq of 8 Tesla is driven by quadrupolar 
fluctuations and that dipolar correlations are short-range with moments parallel to the applied 
magnetic field in the high-field phase. Heat-capacity measurements evidence a phase transition 
into this high-field phase, with an anomaly clearly different from that at low magnetic fields. Our 
experimental data are consistent with a picture where the ground state above Hq has a next-nearest 
neighbour bond-nematic order along the chains with a fluid-like coherence between weakly coupled 
chains. 



PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Kt 

Frustrated quantum ferromagnets are predicted to dis- 
play a new exotic state of matter called a bond-nematic 
state, a highly correlated quantum spin state resembling 
nematic liquid crystals [TH5]. Starting from a magneti- 
cally saturated high-field state, this exotic phase arises 
from Bose-Einstein condensation of two-magnon pairs 
upon lowering the magnetic field. Unlike conventional 
Bose-Einstein condensates, there is no transverse long- 
range dipolar order but (quasi-) long-range quadrupolar- 
nematic order. The corresponding order parameter is 
defined on the bond between two neighboring sites by a 
rank-2 tensor product of the local spin degrees of free- 
dom. This new type of quantum ground state should be 
observable in real materials, since it is predicted to exist 
for a large variety of ID [3H7], 2D [TJ d M M and even 
some 3D spin- 1/2 ferromagnets [10] . which all have in 
common a strong next-nearest-neighbor frustrating ex- 
change. A particularly interesting model is the spin-1/2 
frustrated ferromagnetic Heisenberg (FFH) chain, where 
a quadrupolar-nematic phase is predicted to extend from 
low magnetic fields to the saturation field i? S at, with al- 
gebraic decay of transverse quadrupolar and longitudinal 
dipolar correlations [3HS1 UJ] ■ In view of the strong theo- 
retical activity in this field, clearly identified experimen- 
tal realizations of bond-nematic states are needed. 

The FFH chain LiCuVC>4 is an excellent candidate to 
search for an experimental realization of bond-nematic 
correlations. The excitation spectrum displays pro- 
nounced quantum effects [121 113] and a magnetization 
anomaly near fioHc ~ 40-48 T [T3] corresponds to a 
predicted long-range ordered quadrupolar bond-nematic 
phase p]. Between Hq and an additional magneti- 
zation anomaly at HqHq ~ 8 T [15], the order has 
been supposed to be dipolar long-ranged and amplitude- 



modulated pi)] IT7] . However, the character of the or- 
dered phases above Hq remains a subject of speculation. 

In this Letter we present experimental evidence that 
in LiCuVC>4 the field-induced phase transition at Hq is 
driven by quadrupolar correlations. Based on neutron 
scattering, we show explicitly that the dipolar correla- 
tions, which are long-range (LR) at zero field due to a 
weak interchain coupling, are in fact short-range (SR) 
in all crystallographic directions above Hq down to the 
lowest temperatures (100 mK) and involve only the spin 
components parallel to the magnetic field. The field- 
dependent characteristic wave vector follows precisely 
the theoretical predictions for the quadrupolar-nematic 
phase, namely fcic = (1/2 — m)/p, where m = (S*) is 
the uniform magnetization and p — 2 for the quadrupolar 
bond-nematic state [3]-|6]. A concomitant change in or- 
der parameter (or universality class) is revealed by heat 
capacity measurements. 

The Cu 2+ ions in LiCuVC>4 form spin-1/2 chains 
along the orthorhombic b axis [HI HH] • Nearest-neighbor 
(NN) spins are coupled ferromagnetically (FM) ( Ji < 0) 
and next-nearest neighbors (NNN) antiferromagnetically 
(AF) (J2 > 0) [12 . Small interchain interactions lead 
to 3D LR order below T/v ~ 2.4 K, where the reduced 
ordered moments (mo ~ 0.3/zs) form a circular cycloid 
in the a6-plane with incommensurate propagation vec- 
tor fcic = (0,0.532,0) [H CS3, which corresponds to 
k\c = 0.468 r.l.u. for a ID chain due to the centering 
translation of the Imma space group and fcic = 0.234 in 
units of 2-7r/<iNN, where c?nn = b/2 is the NN distance. 

A first series of unpolarized elastic neutron scattering 
measurements were performed on the triple-axis spec- 
trometers PANDA at FRM II (Munich) and IN14 and 
IN20 at Institut Laue-Langevin (Grenoble), with verti- 
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FIG. 1. (Color online) (a-b) Incommensurate propagation 
vector fcic in units of 2-7r/dNN as a function of applied field 
for H || c (left) and H || a (right) compared to the theoretical 
prediction for the quadrupolar-nematic phase [3HB] , using the 
experimental magnetization m a , c of Ref. [14] . normalized to 
1/2 at the saturation field, (c-d) Specific heat for H||c (left) 
and H || a (right). The distinct anomaly evidences a phase 
transition at all fields and the change in shape near Hq ~8 T 
suggests a change of the universality class. 



cally focusing monochromator and analyzer and incident 
wave vectors 1.5 < fcj < 2.662 A -1 . Second order con- 
tamination was suppressed by a cooled Be-filter or a py- 
rolythic graphite filter. The same single crystal as used 
in Ref. |13j was mounted with either the a- or the c-axis 
along the vertical field of a 15 T split-coil cryomagnet. A 
dilution insert cooled the sample to 100 mK. 

Figure [lja,b) shows the magnetic field dependence of 
the propagation vector of dipolar spin correlations within 
the chain, kic- Below (j,qHq^8 T, fcic is essentially field- 
independent while above Hq it follows the relation kic = 
(1/2 — m a , c (H))/p with p — 2, where m a and m c are the 
magnetization curves measured at 1.3 K |14j . normalized 
to m = 1/2 at the saturation field. Both below and 
above Hq, k\c behaves as predicted theoretically [SHE]. 
Above Hq the slope of kic versus m, i.e. p = 2, proves 
the presence of quadrupolar correlations for both H || c 
[Fig. 0a)] and H || a [Fig. [gb)]. 
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FIG. 2. (Color online) (a-c) Scans of the dipolar correlations 
at T = 0.1 K for the three reciprocal space directions are reso- 
lution limited at zero field and broadened at high fields above 
Hq ~8T. (d-f ) Magnetic field dependence of the inverse dipo- 
lar correlation length k = The onset of short-range order 
at Hq is clearly seen for all three reciprocal space directions. 



The specific heat close to the phase transition from the 
paramagnetic state into the magnetically ordered phases 
was measured employing the relaxation method using a 
Quantum Design calorimeter. Despite an obvious change 
in the peak shape, see Fig. [lTc,d), the anomalies at 9 T 
for H || a and H || c remain remarkably well-defined and 
clearly different from an SR order maximum, thus testi- 
fying a phase transition. The pronounced change of the 
shape of the specific heat peak above Hq indicates a dif- 
ferent universality class compared to lower fields. The 
peak shape above Hq is identical for H || c and H j| a, 
evidencing the same universality class for the two field 
directions. 

Below Hq, the order parameter is dipolar as shown by 
the sharp (resolution-limited) magnetic Bragg peaks in 
our neutron scattering measurements. Above Hq, these 
peaks are clearly broadened for all three crystallographic 
directions, as shown in Fig. |2ja-c). The increase in the 
peak width appears abruptly at Hq, as evidenced in 
Fig. [2jcl— f ) , where the intrinsic peak width k, which is 
inversely proportional to the dipolar correlation length 
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FIG. 3. (Color online) Polarized cross-sections measured at 
T = 70 mK for the magnetic reflections Q = (1, kic, 0) with 
H || c (left panels, a-c) and Q = (0, — kic, 1) with H || a (right 
panels, d-f). 



£, is plotted as a function of applied field for the three 
space directions. The occurrence of a finite dipolar cor- 
relation length £ above Hq clearly establishes that the 
order parameter is no longer dipolar. In order to rule out 
the hypothetical existence of a dipolar LR order with a 
different wave vector, we have estimated the integrated 
intensity involved in the SR dipolar correlations. We 
first note that both fcic an d Hk are ^-independent, since 
scans along (0, k, 1.5) and (0, k, 1) for H || a yield iden- 
tical kic and Kfc. The integrated intensity above Hq 
can then be obtained by explicit integration over £ in 
(0,fci C ,£) [Fig. go)] and k in (0, k, 1) [Fig. gb)], while 
the /i-integration is effectively achieved by the coarse ver- 
tical resolution in the measurements. The estimated in- 
tegrated intensity at 14.5 T is at least as high as at zero 
field, i.e., the dipolar SR correlations at high field in- 
volve at least as much of the magnetic moment as the 
LR order at zero field. We can therefore exclude a con- 
ventional dipolar LR order with another ordering wave 
vector. The disparate correlation lengths of the SR dipo- 
lar correlations along a, b, and c may be related to the 
hierarchy of the exchange interactions in LiCuV04 [T2] • 
To determine which spin components are involved in 
the SR dipolar correlations above Hq, we performed a 
second series of elastic neutron scattering experiments, 
employing polarized neutrons in a vertical magnetic field 
with H || c or H || a. The measurements were performed 
on the cold triple-axis spectrometer IN14 using an asym- 
metric cryomagnet reaching up to 12 T and a dilution 
refrigerator going down to 70 mK. Neutrons of wave vec- 
tor 1.5 A -1 from a PG(002) monochromator were po- 



larized by a supermirror bender and their polarization, 
being parallel or anti-parallel to the applied field, was 
controlled by a spin flipper in the incoming beam and 
analyzed with the (111) reflection of a Heusler crystal. 
After flipping ratio corrections and subtraction of the 
spin-incoherent background, the measured non-spin-flip 
(NSF) and spin-flip (SF) cross-sections are in this set-up 
proportional to dipolar correlations polarized parallel to 
the field (perpendicular to the scattering plane) and per- 
pendicular to the applied field (and also perpendicular to 
the wavevector transfer Q). The results, shown in Fig.|3j 
will now be discussed. 

The magnetic Bragg reflection (l,/cic,0) measured for 
H || c [Fig. |3][a-c)] evolves from purely SF character at 
2 T to purely NSF character at 9 T with a remarkable 
coexistence of the two components in the vicinity of Hq 
at 7.8 T. This shows unambiguously that the dipolar 
moments swap in a first order transition from the ab- 
plane, transverse to the applied magnetic field, to the c 
axis, parallel to the applied field. Similarly, the reflection 
(0, — kic, 1) measured for H || a above the 3.5 T spin-flop 
field [Fig.[3]jd-f)] reveals that the scattering changes from 
SF at 7 T (transverse dipolar, now with moments in the 
be plane) to NSF at 9 and 10 T with moments along the 
field direction a. For both magnetic-field orientations we 
therefore conclude that the incommensurate dipolar cor- 
relations are LR and transverse to the applied field at low 
fields and become SR and longitudinal for fields above 
Hq. The use of neutron scattering allows us to draw 
such conclusions without an a priori model. Since the 
dipolar correlations become SR above Hq, the distinct 
anomaly in C p /T above Hq cannot result from a dipolar 
order parameter. The particular field dependence of the 
characteristic wave vector fcic shown in Fig. [lja,b) im- 
plies that the phase transition is driven by quadrupolar 
correlations. 

We will now compare our experimental results on 
LiCuV04 to theoretical predictions. We start with the 
purely ID 5 = 1/2 FM frustrated Heisenberg (FFH) 
chain with NN FM exchange Ji, AF NNN exchange J 2 , 
and IJ1/J2I < 2.5 [SH3 [H]- At zero and small mag- 
netic field H || z, uniform LR order of the bond vector- 
chirality (S n x S n+ i) z is predicted while transverse and 
longitudinal dipolar correlations are critical with a field- 
independent propagation vector close to k — 1/4, indi- 
cating incipient spiral order E]. The magnetization 
increases in steps of AS Z = 1 [4H3 [21]. These charac- 
teristic features are illustrated using a soliton (fermion) 
picture in Fig. Qa) , where the magnetization is increased 
by adding S = l two-soliton bound states into the ground 
state. Hq is identified as a transition to AS Z = 2 
steps in the magnetization [3l [5l |21] [Fig. Qb)]. This is 
the signature for a two-magnon Bose condensate ground 
state [H [5], with critical transverse bond-quadrupolar 
correlations (Sq S^ n+1 S^ S^ +2m+1 ) [S], including cor- 
relations of the two-magnon bound state propagator, 
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FIG. 4. (Color online) (a) FFH-chain at H — seen as two 
intercalated Heisenberg AF chains (closed and open arrows), 
vertically displaced to illustrate the direction of the NN-bond 
vector chirality (green crosses and red dots). Neel solitons 
(red solid arcs) propagate freely on a Heisenberg AF chain 
and are bound to S — 1 states by Ji in the FFH-chain, result- 
ing in vector-chiral LR order. (b,c) FFH-chain for H > Hq. 
(b) Two-magnon (four-soliton) states provide AS Z = 2 mag- 
netization steps, (c) Propagating individual solitons costs no 
extra Ji-energy, as long as one soliton in the even subchain 
is followed by one in the odd subchain. This non-local order 
is reflected in the AF sequence of the bond pseudo-spin 
(numbers). (d,e) Two FFH-chains coupled by a FM diagonal 
interchain interaction J5. (d) H — 0: J5 binds two two-solitons 
into a four-soliton bound state and thus stabilizes dipolar LR 
order, (e) H > Hq: Propagation of individual solitons costs 
no extra Js-energy, as long as the non-local AF order of the 
bond pseudo-spin is kept. 



(S^S^S^S W ) [MS]. The density fluctuations of the 
two-magnon bound state condensate [U |5j lead to criti- 
cal correlations (S^S^) [HHS] with a propagation vector 
k = 1/4 — to/2, characteristic for the whole quadrupolar- 
nematic two-magnon condensate phase. Transverse dipo- 
lar correlations are SR in this phase [3H5]. These char- 
acteristics can likewise be illustrated in a soliton image 
[Fig. [4jc)] , where propagation of individual solitons de- 
stroys dipolar LR order. The predictions for the perfect 
FFH-chain correspond well to the observation of vector- 
chiral LR order in LiCuV04 |20j with a field-independent 
k\c at low field, and the presence of longitudinal dipo- 
lar correlations with a propagation vector characteristic 
of quadrupolar correlations above Hq, observed in this 
work. 

However, LiCuVC>4 does not correspond exactly to the 
purely ID FFH-chain discussed above, since it has a weak 



easy-plane anisotropy [2U] and weak interchain interac- 
tions [12] . We now consider these two effects in turn. 
The main effect of exchange anisotropy is to shift Hq 
to higher fields [TJ, without any qualitative changes to 
the vector-chiral and quadrupolar phases in the range 
of J1/J2 relevant for LiCuVC>4. However, the large 
value of Hq of 8 T cannot be reconciled with the very 
weak easy-plane anisotropy of only 0.5% [20] . The ef- 
fect of 2D interchain interactions has been studied in the 
high-field isotropic exchange limit [5]. Close to the sat- 
uration field -ff sa t, the transverse quadrupolar correla- 
tions \Sq Si S~[ +l ) are shown to be LR ordered, while 
(apart from the uniform component m) the longitudinal 
and transverse dipolar correlations are both SR. Below a 
critical field H c (close to H S3jt ), this LR ordered bond- 
quadrupolar phase becomes unstable. In LiCuVCU, Hq is 
found at 40-48 T from magnetization measurements [2] . 
Surprisingly, LiCuVC-4 shows the essential features of this 
high-field bond-quadrupolar phase already just above 
Hq, namely SR dipolar correlations, quadrupolar corre- 
lations, and a phase transition evidenced by the specific 
heat. Our work does not give support for a hypothetical 
LR longitudinal dipolar (amplitude-modulated) phase in 
the lower-field region of the quadrupolar-nematic phase 

dHzi ng na E2] - 

We now argue that the peculiar type of order observed 
in LiCuV04 is related to the diagonal nature of the dom- 
inant interchain exchanges: the FM J5 acts diagonal in 
the ab plane, while the much smaller AF J§ acts along 
the main body diagonal [12]. Within the soliton pic- 
ture, at zero magnetic field, J5 and Jq favor bound four- 
soliton states, and thus spiral dipolar LR order, which 
is only little perturbed by bound four-soliton propaga- 
tion [Fig. Qd)]. This completely agrees with the ob- 
served dipolar LR order in LiCuV04 at H — 0. More- 
over, inelastic neutron scattering evidences bound two- 
spinon (two-soliton) excitations and additional strong 
four-spinon excitations 13], which may find their expla- 
nation in a ground state that already contains bound 
four-solitons. Above Hq, in contrast to the zero-field 
case, neither J 5 nor J 6 bind the four solitons [Fig. Qe)], 
and the propagation of individual solitons destroys lon- 
gitudinal dipolar LR order in all reciprocal space direc- 
tions, as we observe experimentally in LiCuV04. 

In the soliton picture of the FFH-chain above Hq, the 
NNN bond pseudo-spins, a z n = l(-l) n (4S*S* +2 + 1) 
(whose values are reported in Fig. |4j), form a non-local 
LR AF string order with finite values of the correla- 
tion function G(n) = « o expjm Emt™„° ? <J<+„ D ) for 
n — > 00 [23] . resembling the string order in valence-bond 
solid spin-1 chains and the Haldane ground state. This 
means that for finite FM J\, a propagating pseudo-spin 
afj = +1 cannot go past neighboring of„ = — 1 and vice 
versa. Since only the direction of the bond-pseudo spin 
is ordered, but not its position, the AF string order of 
the bond-pseudo spin is nematic in a sense very close 



5 



to liquid-crystalline phases. Finite FM J5 or AF Jq not 
only favor the AF string order along the FFH-chains, 
but also coherence between chains such that a pseudo- 
spin er* = +1 cannot propagate past er^ = — 1 states in 
neighboring chains. This leads to SR dipolar correlations 
in the a6-plane and between planes. Since the LR string 
order breaks at least one discrete Z2 symmetry, with a 
certain 2D or 3D coherence, it could lead to a finite- 
temperature phase transition. We speculate that the dis- 
tinct anomaly in the specific heat above Hq evidences a 
transition to such a 2D or 3D correlated low-temperature 
phase with AF alignment of the bond pseudo-spins but 
no positional order. Further theoretical investigations of 
the order parameter in the quadrupolar-nematic phases 
of coupled FFH-chains are clearly needed to fully explain 
our experimental results. 
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